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ABSTRACT 


The  performances  of  the  conventional  coherent  beam  envelope 
detector  and  the  incoherent  beam  envelope  detector  were  studied  and  compared 
using  short-period  Norwegian  Seismic  Array  (NORSAR)  data  from  regional 
earthquakes  and  presumed  explosions.  The  detectors  were  evaluated  in  two 
frequency  bands  in  terms  of  detection  performance,  false  alarm  characteris- 
tic, and  the  resultant  operating  characteristic.  On  the  basis  of  the  derived 
operating  characteristics  from  data,  the  detectors  with  specific  operating  pass- 
band  ai  e ranked  in  the  capability-decreasing  order;  incoherent  detector 
(1. 5-2.  5 Hz),  coherent  detector  (1.  5-2.  5 Hz),  incoherent  detector  (3.  0-4.  0 Hz), 
and  coherent  detector  (3.  0-4.0  Hz).  A detector  using  the  envelope  of  the  array 
beam  computed  with  the  Hilbert  transform  was  also  examined. 


Neither  the  Advanced  Research  Projects  Agency  nor  the  Air  Force 
Technical  Applications  Center  will  be  responsible  for  information  contained 
herein  which  has  been  supplied  by  other  organizations  or  contractors,  and  this 
document  is  subject  to  later  revision  as  may  be  necessary.  The  views  and  con- 
clusions presented  are  those  of  the  authors  and  should  net  be  interpreted  as 
necessarily  representing  the  official  policies,  either  expressed  or  implied,  of 
the  Advanced  Research  Projects  Agency,  the  Air  Force  Technical  Applications 
Center,  or  the  US  Government. 
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SECTION  I 
INTRODUCTION 


I 

I 
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Tliis  report  presents  the  rei'ults  of  a study  of  detection  on  the 
envelope  of  "ihort-period  array  beams  from  the  Norwegian  Seismic  Array 
(NORSAR)  which  had  been  formed  through  either  incoherent  summation  or  conven- 
tional coherent  summation  techniques.  Coherent  beamforming  of  a large  short-period 
array  does  not  yield  the  theoretically  predicted  signal-to -noise  ratio  improve- 
ment due  to  imperfect  signal  similarity  across  the  array.  This  signal  dissimi- 
larity is  frequency  dependent  with  the  moat  pronounced  drop  in  signal  coherence 
occurring  at  frequencies  above  ZHz,  Thus,  high  frequency  signals  suffer  the 
largest  losses  in  signal-to-noise  ratio  (SNR),  At  NORSAR,  this  type  of  signal 
originates  from  near-regional  events  and  Eurasian  presumed  explosions,  par- 
ticularly those  from  Western  Russia  (Ringdal  et  al.  , 197Z). 

Signal  loss  also  may  occur  even  with  well -equalized  and  similar 
signal  waveforms  if  the  regional  corrections  are  not  known  or  if  the  signal  ar- 
rives slightly  off-azimuth  from  the  direction  of  a preformed  beam.  Whenever 
signal  dissimilarity  is  severe  enough  to  cause  appreciable  loss  in  signal-to- 
noise  ratio  during  conventional  coherent  beam  forming  or  when  beaming  para- 
meters are  in  slight  error,  incoherent  beamforming  using  the  sum  of  the  rec- 
tified sensor  or  subarray  output  may  allow  partial  recovery  of  the  signal  loss. 

The  objective  of  this  study  was  to  investigate  systematically 
the  performance  of  the  envelope  detector  in  detecting  underground  explosions 
and  near-regional  earthquakes.  The  following  steps  were  undertaken  to  ac- 
complish this  objective: 


• Measuring  the  false  alarm  probability  of  coherent  and  incoherent 
envelope  detectors. 

• Computing  the  detection  probability  of  coherent  and  incoherent 
envelope  detectors  for  presumed  underground  explosions  and 
near-regional  earthquakes. 

• Comparing  the  operating  characteristics  of  the  detectors. 

Section  II  describes  the  detection  algorithms  and  the  mathe- 
matical basis  of  the  false  alarm  and  detection  probabilities.  Section  III  de- 
scribes the  data  base.  Section  IV  discusses  the  results  of  this  study.  Finally, 
a concluding  summary  is  made  in  the  last  section. 
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SECTION  II 


THE  dete::tion  algorithms  and  the  fundamentals 

OF  FALSE  ALARM  AND  DETECTION  PROBABILITIES 


A.  INTRODUCTION 

This  section  describes  the  mathematical  operations  of  the  detec- 
tion system  and  the  general  concept  of  false  alarm  and  detection  probabilities 
used  in  evaluating  detection  performance.  The  coherent  and  incoherent  array 
beamformers  are  described  in  Subsection  B,  followed  by  a description  of  en- 
velope forming  in  Subsection  C.  Subsection  D describes  the  detector  output. 
The  last  three  subsections  introduce  the  fundamentals  of  false  alarm  and  detec- 
tion probabilities,  and  detector  operating  characteristics,  respectively. 


B.  COHERENT  AND  INCOHERENT  BEAMFORMING 

Let  S^(t)  be  the  time  trace  of  subarray  i . The  conventional 
coherent  beam  is  defined  as 


ABc(t) 


(II-l) 


where  AB  (t)  is  the  coherent  array  beam,  t.  is  the  travel-time  delay  of  the 
C-*  i 

signal  waveform  for  subarray  i and  N is  the  number  of  channels. 


ABj(t)  . 


Following  Ringdal  et  al.  , (1972),  the  incoherent  array  beam, 
was  formed  by  summing  the  absolute  values  of  the  subarray  beams: 


ABj(t) 


(II-2) 


II- 1 


where  the  symbols  used  here  have  the  same  meaning  as  in  equation  (II- 1). 


In  the  work  described  here,  each  channel,  (i.  e.  , subarray 
beam  S.(t)  ),  was  prefiltered  with  a zero-phase  bandpass  filter  before  forming 
the  array  beam.  To  compare  detector  performance  in  different  frequency 
ranges,  two  passbands,  1.  5 -2.  5 Hz  and  3.  0-4.  0 Hz,  were  used.  Both  sides 
of  each  passband  were  given  a 0.  7 Hz-wide  cosine  taper.  Figure  II- 1 shows 
the  filter  amplitude  response  for  the  two  passbands. 


Both  adjusted-delay  and  diversity-stack  techniques  were  used 
for  array  beamforming.  The  travel  time  delay  t.  for  each  channel  appro- 
priate for  each  signal  was  computed  either  from  the  cross  correlation  function 
with  a reference  channel  or  from  manual  input.  For  adjusted-delay  beamforming, 
each  channel  is  equally  weighted;  for  diversity-stack  beamforming,  each  chan- 
nel was  weighted  by  the  quantity  W(i)  which  was  computed  according  to  the 
formula  (Texas  Instruments,  1971): 


W(i) 


(i)  - P„(i) 
npii) 

n 


(II-3) 


where  P (i)  and  P (i)  are  signal  and  noise  powers  for  channel  i , respec- 
s n ~ 

lively. 


C. 


THE  ENVELOPES 


1.  Short-Term- Averages  As  Envelopes 


The  envelopes  of  the  coherent  and  incoherent  array  beams  AB  (t) 

C 

and  ABj(t)  defined  by  equations  (II-l)  and(lI-2)  were  formed  by  taking  a running 
average  over  a specified  gate  length: 


(II-4) 


/=0 


where  STA  refers  to  short-term-averages  and  is  used  to  represent  the  en- 
velope amplitude  of  signals,  At  is  the  sample  interval  and  L+1  is  the  num- 
ber of  points  averaged.  For  the  short-period  NORSAR  data  processed,  At 
was  0.  1 seconds  and  L+1  was  set  to  15,  equivalent  to  1.5  seconds,  which  is 
the  timegate  used  in  the  NORSAR  on-line  detection  system. 

The  STA  envelope  is  fast  and  straightforward  to  compute;  how- 
ever, strictly  speaking,  it  is  only  an  approximation  of  the  envelope  of  a signal 
wave  train.  The  envelope  of  a time  series  is  defined  by  the  magnitude  of  its 
analytic  function  computed  from  its  Hilbert  transform. 


2.  Amplitudes  of  Analytical  Signal  As  Envelopes 

Since  the  short-term-average  represents  only  the  approximate 
shape  of  a wave  train  envelope,  the  maximum  amplitude  of  a signal  is  not  ac- 
curately measured  in  the  STA  representation.  Because  the  true  envelope  am 
plitude  may  allow  better  bodywave  magnitude  measurements  and  may  also  in- 
crease the  detectability  of  event  signals,  an  attempt  was  made  to  form  the 
exact  envelope  of  a wave  train  by  Hilbert  transformation. 

The  Hilbert  transform  of  a time  series  x(t)  is  the  time  series 
x(t)  defined  by  the  convolution  integral  (Cizek,  1970): 


x{t)  = 


00 

-L  f 

IT  J r - 


dr  * x(t)  . 

7Tt 


(II-5) 


Denoting  X(f)  as  the  Fourier  transform  of  x(t)  and  H(f)  as 


that  of  ( ),  we  have 

7Tt 


where 


F |x(t)|  = X(f)  . H(f) 

( -j  if  f < 0 

H(f)  = Jo  if  f=0 

* +j  If  f > 0 


(II-6) 


and 


1 


I i 

I « 


ii 

n 

n 

n 

II 

n 

!i 

;i 


To  compute  equation  (II-5)  in  discrete  form,  let  x.  , i = 0,  1, 


discrete  Fourier  transform  (DFT)  of  this  sequence  is  a sequence  X , k=0, 

k 


M-1,  denote  the  sequence  of  M complex  finite  values  of  x(t).  Then  the 
screte  Fourier  transfc 
1 , . . . , M-  1 , defined  by 

- _L  V ^ ‘ j il<(2  7r/M) 

M / J i 


M-1 


X. 


(II-7) 


i=0 


For  an  inverse  DFT,  the  sequence  x.  is  defined  by 


M-1 


E 


+ j i k (Ztt/  M) 


(11-8) 


k=0 


Equation  (Il-‘i)  holds  true  for  the  Fourier  transform  of  a Hilbert 

transform.  Therefore,  the  discrete  Hilbert  transform  can  be  performed  in  the 

frequency  domain  through  Fourier  transformation.  In  doing  so,  Jet  G be  the 

k 

sequence  of  equation  (11-6): 


(11-9) 


whe  re 


i -j 

if 

k = 

H = 
k 

1 " 

if 

k = 

0.  M/2 
M 

+J 

if 

k = 

T ^ 

The  discrete  sequence  of  the 

p roc 

«*ss 

x(l)  is 

M 

2 


the  sequence  G in  equation  (H-9). 


Define  the  coin])l<-x  process  /,(t)  by  taking 


) — 

t ( 

[ J. 

I I 
t 


/.(t)  = x(t)  + Jx(t) 


(H-10) 


li-h 


This  process  is  known  as  the  analytic  signal  assj^.ciated  with  x(t)  . Finally, 
the  envelope  of  x(t)  is  the  amplitude  of  z(t)  (Bracewell,  1965)  which  is 


E(t)  • [x^t)  + . 


(II- U) 


For  coherent  array  beamforming,  x(t)  in  equation  (II-ll)  is 
obtained  from  AB^(t)  in  equation  ^11- 1),  The  incoherent  array  envelope  Ej(t) 
was  formed  by  averaging  e suba^ray  envelopes: 


(11-12) 


where 


is  the  envelope  of  subarray  i , 
is  the  wave  travel- time  delay,  and 
N is  the  number  of  subarrays. 


D.  THE  DETECTOR  OUTPUT 

In  the  decision-making  process  of  a detection  system,  the  esti- 
mated signal-to-noise  ratio  (SNR)  commonly  is  used  as  the  detection  statistic. 

A detection  is  declared  when  the  SNR  exceeds  a preselected  decision  threshold. 
Here,  signal  levels  were  measured  by  the  quantities  STA  and  E.  For  opera- 
tional detectors,  the  conventional  computation  of  noise  RMS  values  is  somewhat 
tedious  and  is  difficult  to  realize  in  practice.  Thus  the  Long-Term-Averages 
(LTA)  was  introduced  to  represent  the  noise  level.  The  LTA  is  computed  in  the 
same  manner  as  in  equation  (II-4),  but  with  an  integration  time  of  30  seconds. 
Because  the  LTA  is  averaged  over  a much  longer  time  window  than  the  STA,  it 
is  a good  approximation  of  the  noise  mean  value.  The  signal-to-noise  ratio 
(expressed  in  dB)  for  the  detection  system  used  here  was  thereafter  denoted  by 


II-6 


(11-13) 
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SNR(t)  = 20  X logjQ  (STA(t)/ LTA(t)) 

STA  -nust  be  replaced  with  E(t)  for  the  envelope  computed  from  the  Hilbert 
transform. 

The  process  in  equation  (11-13)  was  the  actual  operational  mode 
in  event  detection.  However,  the  running  LTA(t)  was  contaminated  with  signal 
energy  when  maximum  SNR(t)  occurred.  To  obtain  a more  meaningful  estimate 
of  the  signal-to-noise  ratio,  the  LTA  was  frozen  prior  to  the  signal  arrival,  the 
time  of  which  was  determined  from  the  event  location  and  origin  time.  The 
maximum  amplitude  of  the  envelope  in  a signal  gate  was  used  for  the  computa- 
tion of  signal-to-noifce  ratio: 

SNR  = 20x  logjQ((STA)^^^/LTA)  (11-14) 

where  STA  is  replaced  with  E if  the  analytic  signal  is  used,  and  where  LTA 
is  the  average  of  LTA  in  the  noise  gate.  Both  processes  of  equations  (11-13) 
and  (11-14)  were  performed. 


E.  FALSE  ALARM  PROBABILITY 

When  a detector  output  exceeds  the  decision  threshold  while 
there  is  no  signal  present,  it  is  a false  alarm  known  as  Type  I Error  in  the 
detection  problem.  False  alarm  probability  can  be  measured  by  forming  the 
histogram  of  detector  output  in  noise.  This  normalized  histogram  is  essen- 
tially the  false  alarm  probability  density  function  for  the  detection  system.  If 
the  noise  amplitudes  are  log  normally  distributed,  the  probability  density  func- 
tion follows  a Gaussian  distribution: 


dF 

dS 


Pf(S)  = 


(11-15) 


• » 
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where 


S is  the  detector  output  for  noise, 

is  the  mean  value,  and 
n 

cr  is  the  standard  deviation, 
n 

Thus,  dF/dS  or  P^(S)  is  the  false  alarm  probability  density  function. 

False  alarm  probability  for  a given  decision  threshold  S*  is 
obtained  by  integrating  the  density  function  from  S*  to  infinity,  or,  mathema 
tically 

00 

P£(J?)dS.  (11-16) 

S* 

This  cumulative  distribution  function  can  be  used  to  determine  the  probable 
number  of  false  alarms  registered  by  the  detector  for  a given  time  period. 


F.  DETECTION  PROBABILITY 

When  the  decector  fails  to  detect  a signal  for  various  reasons 
such  as  signal  dissimilarity,  erroneous  look  direction,  or  weak  events,  the 
failure  is  called  a Type  11  Error.  Assuming  the  signal  amplitude  to  be  a ran- 
dom variable  the  probability  of  detecting  a signal  at  a given  level  is  called  the 
detection  probability  for  that  level. 

If  the  signal  amplitude  is  log  normally  distributed,  the  probabil 
ity  density  function  of  the  output  level  follows  a Gaussian  distribution: 


dD 


1 


dS  = Pd'"’  = 

V Ztt  cr 


exp 


as  its  cumulative  function  from  to  infinity 


D = 


e 

f 


P AS)  dS 


(11-17) 


(11-18) 


I 

I 

I 


is  thus  the  probability  of  detection  at  the  decision  threshold  S*, 
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Experimentally,  the  detection  probability  density  function  is  ob- 
tained by  forming  the  histogram  for  the  detector  output  from  events  and  the  de- 
tection probability  by  integrating  the  density  function  from  a given  level  to  in- 
finity. 

G.  DETECTOR  OPERATING  CHARACTERISTIC 

From  the  false  alarm  probability,  F(S),  and  the  detection  proba- 
bility, D(S)  of  a given  detector,  the  common  factor  S can  be  eliminated  to  form 
a new  function  S(F,  D)  which  is  known  as  the  Detector  Operating  Characteristic 
or  Detector  Operating  Curve. 

The  operating  curve  is  a display  of  detection  probability  versus 
the  false  alarm  probability  for  a given  detector -output  level  (SNR).  For  vari- 
ous detectors,  this  display  enables  us  to  evaluate  the  relative  detector  perfor- 
mance by  comparing  the  detection  probabilities  for  a given  false  alarm  rate  or 
the  false  alarm  rates  for  a given  detection  probability.  The  detector  curves  de- 
rived from  the  data  will  be  shown  in  Subsection  IV-E. 
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SECTION  III 


DATA  BASE 


The  data  used  in  this  report  were  NORSAR  short-period  subar- 
ray beams  of  near-regional  earthquakes  and  presumed  underground  explosions 
which  had  been  formed  during  a previous  short-period  NORSAR  array  evalua- 
tion program  (Barnard  and  Whitelaw,  1972;  Ringdal  and  Whitelaw,  1973).  The 
available  data  for  this  purpose  were  very  much  limited.  In  order  to  obtain  a 
more  reliable  statistical  analysis,  all  earthquakes  with  epicentral  distance  less 
than  twenty-seven  degrees  and  all  presumed  underground  explosions  which  were 
available  on  existing  subarray  beam  tapes  were  processed.  In  all,  there  were 
55  earthquakes  and  36  presumed  explosions  which  were  used.  Table  III-l  lists 
the  bodywave  magnitude,  location,  depth,  epicentral  delta  and  source  time  of 
these  events.  In  the  table,  the  first  three  characters  of  event  identification  indi- 
cate the  seismic  region  of  ?.n  event.  The  middle  three-digit  number  indicates 
the  Julian  day  and  the  last  two-digit  number  indicates  the  hour  of  event  time, 
Theother  symbols  or  characters  are  only  for  convenience  in  processing.  The 
distribution  of  the  earthquake  magnitude  is  relatively  narrow  and  is  not  repre- 
sentative of  the  general  seismicity.  The  magnitude  distribution  of  the  presumed 
explosion  is  even  more  narrowly  restricted  and  ersentially  does  not  overlap  the 

earthquake  distribution.  Therefore,  the  m,  -values  distribution  for  the  ensemble 

b 

is  not  statistically  well-behaved. 

The  detector  output  signal-to-noise  ratios  for  both  adjusted  delay 

and  diversity-stack  beams  are  tabulated  in  Table  III-E.  Both  the  updated  LTA, 

(STA/LTA)  and  the  averaged  LTA,(STA)  /LTA  are  shown  in  two  passbands 
max  max 

The  false  alarm  probability  was  calculated  for  ten  noise  samples 
each  320  seconds  long.  Tabid  III-3  describes  these  samples. 
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TABLE  III-l 

EVENT  INFORMATION 
(PAGE  1 OF  2) 


Event  I.  D. 

"'b 

Lat. 

Long. 

Depth  A 

Day 

Time 

NZH/24 1/C5N 

6.3 

73. 3N 

55.  IE 

C.C 

20.5 

08/28/72 

05. 59.57 

NZH/270/06N 

6.4 

73. 4N 

55.  IE 

0.0 

20.5 

09/27/71 

05.  59.'^'^ 

SWR/277/08N 

5.8 

46. 8N 

45. OE 

0.0 

24.1 

10/08/72 

08. 59.58 

URA/082/06N 

5.6 

61 . 3N 

56.  5E 

0.0 

21.6 

03/23/71 

06. 59.56 

KAZ/356/0fiN 

6.0 

47. 9N 

48. 2E 

o.c 

24.8 

12/22/71 

06.59.56 

URA/191/16N 

5.3 

64. 2N 

55.  2E 

o.c 

20.3 

07/10/71 

16. 59.59 

JAP/063/07N 

5.2 

40. 7N 

143. 5E 

0.0 

71.4 

03/04/71 

07.48.39 

KUR/142/12N 

3.  8 

49. ON 

154. OE 

0.0 

66.2 

05/22/71 

12.59.53 

CAN/310/22N 

6.8 

51 .5N 

179. IE 

2.0 

67.3 

1 1/06/71 

22.00.00 

TAD/188/04N 

4.5 

38. 6N 

73.  IE 

0.0 

43.8 

07/C7/71 

03.52.53 

IRA/102/19N1 

6.0 

28. 3N 

55.  6E 

0.0 

44.1 

04/12/71 

19.03. 26 

KDR/077/03N 

5.1 

49. 2N 

156. 3E 

o.c 

66.5 

C3/18/71 

03. 13. 15 

KH1/073/12N 

5.3 

53. 9N 

160. 5E 

0.0 

62.8 

03/14/71 

12.15.14 

NEV/230/1 4N 

5.4 

37.  IN 

1 16.CW 

0.0 

72.0 

08/18/71 

14.9C.00 

KAZ/070/04N 

5.5 

49. 8N 

78. 2E 

o.c 

38.0 

03/10/72 

04.56.57 

EKZ/363/04N 

4.5 

50. ON 

78. OE 

o.c 

37.8 

12/28/72 

C4. 26.58 

EK2/345/04M 

6.0 

50.  IN 

70. 8E 

o.c 

38.1 

12/10/72 

C4.27.08 

EK1/345/04N 

5.7 

49. 8N 

78.  IE 

0.0 

38.0 

12/1C/72 

04.26.58 

EKZ/246/08N 

5.1 

50. ON 

77. 7E 

o.c 

37.6 

09/02/72 

C3.56.58 

EKZ/239/03N 

5.5 

50. ON 

77. 8E 

o.c 

37.7 

08/26/72 

03. 46.57 

EKZ/229/03N 

5.2 

49. 8N 

78.  IE 

C.C 

38.0 

08/16/72 

03.16.57 

EKZ/188/0 IN 

4.4 

49. 7N 

78. OE 

o.c 

38.0 

07/C6/72 

01.02. 58 

EKZ/041/0SN 

5.5 

50. ON 

78. 9E 

0.0 

38.2 

C2/10/72 

05.02.  ^^7 

EKZ/3f>4/06N 

5.8 

49. 7N 

78.  IE 

0.0 

38. C 

12/30/71 

06. 20.58 

EKZ/349/07N 

4.9 

50. ON 

77. 9E 

o.c 

37.7 

12/15/71 

07.52.59 

EKZ/333/06N 

5.5 

49. 8N 

78.  IE 

o.c 

38. C 

11/29/71 

06.02. 57 

EKZ/294/06N 

5.6 

50. ON 

77. 6E 

0.0 

37.6 

10/21/71 

06.02.57 

EKZ/159/01N 

5.5 

49. 8N 

78. 2E 

0.0 

38.0 

C6/C7/72 

01 . 27.57 

KAZ/18 1/04N 

5.4 

50. ON 

79.  IE 

0.0 

38.0 

06/30/71 

03.56.57 

KAZ/282/06N 

5.4 

50 . 0 N 

77. 7E 

o.c 

37.6 

10/09/71 

06.02.57 

KAZ/157/04N 

5.5 

50. ON 

77. 8E 

0.0 

37.7 

06/06/71 

C4.02.  '^7 

KAZ/145/04N 

5.2 

49. 8N 

78. 2E 

0.0 

38.0 

05/25/71 

C4. 02. 58 

KAZ/1 15/C3N 

5.  9 

49. 8N 

78.  IE 

0.0 

38.0 

04/25/71 

0 3.32.  ■‘8 

KAZ/170/04N 

5.5 

50.  ON 

77. 7F 

0.0 

37.6 

06/19/71 

'■4.0  3.58 

KAZ/081/04N 

5.8 

49. 7N 

78. 2E 

0.0 

38.1 

03/22/71 

04.32.58 

HRS/277/10N 

5.  1 

61. 6N 

47.  IE 

13.0 

17.3 

10/04/71 

10.00.02 

ADS/005/04N 

4.0 

47. 8N 

16.  2E 

11.0 

13.4 

0 1/C5/72 

C4, 57.41 

MRS/248/07N 

4.6 

67. 7N 

33. 4E 

7.0 

1 1.9 

09/04/72 

07.09.04 

WRS/191/07N 

4.6 

52. ON 

31. OF 

3 3.0 

14.  1 

07/09/72 

07. CO. OP 

Ans/189/09N 

4.6 

48. 3N 

14.  5E 

3 3.  C 

12.7 

05/17/72 

9.  2.48 

GPE/200/1 3N 

4.0 

41 . 6N 

23. 8E 

3 3.0 

20.  8 

07/18/72 

13.45.48 

YUG/052/23N 

4.0 

4 1 .ON 

22. 3 E 

’3.0 

21  .0 

02/21/72 

2 3.02. 

CAU/C79/03N 

3.0 

42. 7N 

38.  1’" 

33.0 

24. 4 

03/19/72 

03.34.31 

8IIL/068/22N 

3.5 

40. 8N 

22. 8E 

33.0 

21.3 

C 3/0  8/7  2 

22.04.02 

TUR/21 2/19N 

3.6 

4 1 . 0 N 

27. OE 

33.0 

22.  2 

07/30/72 

19. 46. 24 
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TABLE  III- I 


EVENT  INFORMATION 
(PAGE  2 OF  2) 
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Event  I.  D. 

Lat. 

Long. 

Depth 

A 

Day 

1 

T ime  ] 

AEG/212/01N 

4.4 

39. 9N 

24. 2F, 

33.0 

22.5 

07/30/72 

01.30.0Q 

1 ^ GRE/207/C1N 

4.5 

38. 7N 

21. 4e 

45.0 

23.  1 

07/25/72 

01 . 56.07 

1 

CRT/204/OSN 

4.  1 

47. ON 

36. OE 

33.0 

20. C 

07/22/72 

05.10.47 

1 - GRE/157/10N 

4.2 

37.  ON 

21. 4E 

69.0 

24.0 

06/05/72 

10. 44. 5Q 

1 

GRE/153/1 3N 

4.  1 

39. ON 

24. OE 

33.0 

23.3 

96/01/72 

13.44.11 

1 *•  WRS/202/15N 

3.  8 

52. ON 

54. OE 

33. C 

24.9 

07/21/71 

15.  49.27 

1 . 

SWR/C76/02N 

4.  1 

46. ON 

45.  4E 

33.0 

24.9 

03/17/71 

C2.  1 5. 46 

1 . 

WPS/295/05N 

5.3 

51. 6N 

54. 5E 

6.0 

25.4 

10/22/71 

05.00.C” 

1 ••  GRE/07U/15N 

4.8 

37. 2N 

24. OF 

32.0 

25.0 

03/15/71 

15.23.18 

TUR/0f7/23N 

4.7 

37. 5N 

29. 7E 

21.0 

26.2 

03/08/7 1 

22.44.47 

• 

GRE/1C9/02N 

5.  1 

39. ON 

20.  5E 

16.0 

22.6 

04/19/71 

02.43.52 

i 

• 

; TUP/278/18N 

4.5 

39. ON 

29. 8E 

9.0 

24.8 

10/05/71 

18.53.06 

TUR/276/C7N 

4.7 

38. 9N 

29. 9E 

23.0 

24.9 

10/03/71 

07. 44 . 27 

i| 

, T0R/161/09N 

4.9 

39.  IN 

29. 6E 

3 3.  0 

24.7 

06/10/7 1 

0 9.31.6'^  ' 

TUR/126/0UN 

4.6 

39. ON 

29. 7E 

23.0 

24.  P 

05/96/7 1 

04.24.34  , 

* BLS/210/19N 

4.5 

39. 6N 

30. 4E 

33.  0 

24.4 

C7/29/71 

1 9. 4C  . 1C 

GRE/187/18N 

4.4 

37. ON 

20. OE 

33.0 

24.5 

07/06/72 

18.05.09 

j WRS/179/12N 

3.5 

51. ON 

47. OE 

33.0 

22.  1 

06/27/72 

12.20.36  || 

* GRE/175/07N 

3.4 

37. ON 

21. OF 

33.0 

24.7 

06/23/72 

07.18.14 

TUR/175/04N 

3.7 

41 . ON 

30. OE 

3 3.0 

2 3.0 

06/23/72 

04.25.27 

i 

j TUR/173/05N 

4.  1 

4C . 2N 

30. OF 

33.0 

23.8 

C 6/2 1/7 2 

5.  6.17 

i 

1 TUR/167/14N 

3.3 

38.  ON 

28. OF 

3 3.  0 

25.2 

C6/15/72 

14. 19.02 

GRE/167/00N 

4.9 

38. 3N 

22. 2E 

26.0 

2 3.6 

06/15/72 

C . 33. 24 

• 

ITA/018/23N 

4.  1 

t»4. 2N 

8. 2E 

25.  0 

16.7 

01/18/72 

23. 26.12 

p 

1 

i YUG/067/05N 

2.7 

43.  ON 

21. OE 

3 3.0 

18.8 

03/07/72 

05.21.21 

^ YUG/063/21N 

4.9 

44. 7N 

18. 4E 

33.0 

16.7 

.0  3/C  3/7  2 

21.26.51 

ITA/035/02N 

4.8 

43. 8N 

13.  3E 

2 5.  C 

17.1 

C2/04/72 

02. 4 2.  19 

■ TTA/039/12N 

4.6 

43. 8N 

13.  3E 

33.0 

17.  1 

02/08/72 

12.19. 15 

1 ADR/037/0  IN 

4.9 

44  .ON 

13.  2F 

33.  0 

16.9 

02/06/72 

01.34.22  r 

1 

ITA/036/07N 

4.7 

43. 9N 

13.  3E 

33.0 

17.9 

02/05/72 

07. C8. 13  » 

1 •’  ITA/036/05N 

4.6 

43. 7N 

13.  5E 

33.0 

17.2 

02/C5/72 

05.05.51  1 

1 ..  ITA/C36/03N 

4.4 

43. 2N 

13.  7E 

33.  0 

17.7 

02/C5/72 

03.49.45 

f 

ITA/036/01N 

4.8 

43. 8N 

1 3.  3E 

33.  0 

17.1 

02/05/72 

01.26.23  f 

1 ..  ITA/035/19N 

4.8 

ti3.8N 

13.  3E 

33.^ 

17.1 

02/94/72 

19.02.56 

1 

ITA/035/17N 

4.4 

43. 8N 

13.  3E 

2 3.0 

17.  1 

02/04/72 

17.19.52  i; 

1 •"  ITA/035/C9N 

4.4 

43 . 9N 

13.  2E 

23.  ^ 

17.0 

02/04/72 

09. 18.32  ; 

1 ITA/035/04N 

4.8 

43. 9N 

13.  2E 

33.  0 

17.  0 

02/04/72 

04.40.50  r 

k 

f 

* 

SRR/205/1 IN 

3.5 

48.  ON 

28. OE 

3 3.0 

16.2 

97/24/71 

11.11.42  1 

» YUG/100/02N 

4.6 

42. 5N 

20.  IF 

2 1.0 

19.  2 

C4/10/71 

C2.58.C6  i 

WFS/262/1 IN 

4.5 

57. PN 

41.  IE 

3 3.0 

15.6 

99/19/71 

11.C-9.fO  1 

- ROM/251/04N 

3.4 

45.  BN 

27. OF 

149.  C 

17.8 

C9/63/71 

9 4 . 1 0 . 1 8 

j 

1,  YUG/180/01N 

4.  9 

43. ON 

20. 5E 

33.0 

18.  7 

06/28/72 

01.43.56 

YUG/177/04N 

4.4 

44  . ON 

15.  RE 

33.  0 

17.1 

06/25/72 

04.59.19 

m 

P YUG/176/C7N 

5.3 

43. 7N 

16. 9E 

3 '.  C 

17.5 

06/24/72 

07.17.56 

TTA/173/1 5N 

4.4 

43. 8N 

13.  3f 

4.  0 

17.1 

96/21/72 

15.C6.53  ' 

m 

I 

1 ■ 

^ ITA/166/18N 

[ 

4.Q 

43. 7N 

1 3.  «E 

1 M.  0 
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TABLE  III-2 


DETECTOR  OUTPUT  (dB)  (1.5-2.  5 Hz  BAND) 
(PAGE  1 OF  4) 


Coh.  Incoh.  Coh.  Incoh. 


Event  I.  D. 

A.D. 

D.  S. 

A.E 

>.  D.S. 

A.D. 

, D.S. 

A.  D 

. D.S 

KOP/142/1 2N 

10.2 

8.7 

1.6 

1.9 

12.4 

10.7 

1.8 

2.  3 

KAZ/C70/04N 

24.3 

25.2 

23.0 

23.9 

47.9 

49.9 

37.4 

40.  1 

WRS/277/10N 

22.  7 

22.9 

20.  3 

20.8 

40.  3 

41.4 

33. 6 

34.3 

URA/0P2/06N 

25.3 

25.4 

24.8 

24.8 

63.0 

63.2 

5C.  5 

51.0 

JAP/063/07N 

23.  1 

23.  4 

17.9 

19.0 

35.6 

36.5 

22.  9 

25.0 

EKZ/246/03N 

22.9 

23.4 

18. C 

19.3 

35.  1 

37.0 

22.  9 

25.5 

EKZ/239/03N 

25.  2 

25.  2 

23.  4 

24.  0 

49.8 

51.5 

37.8 

40.3 

EK2/345/04N 

24.4 

25.2 

23.0 

24. 1 

53.0 

54.6 

40.9 

44.6 

EK  1/345/04N 

25.  1 

25.5 

23.9 

24.  5 

53.8 

55.  1 

41.6 

44.7 

SWR/277/08N 

24.8 

24. 8 

23.8 

23.  9 

56.5 

57.0 

45.5 

46.  1 

NZH/241/05N 

25.5 

25.  5 

24.8 

24.  9 

62.2 

62.9 

52.3 

52.7 

EKZ/363/C4  N 

18.8 

19.7 

1 1.C 

12.8 

23.5 

25.  1 

12.  6 

14.8 

EKZ/229/03N 

24.2 

24.8 

22.  1 

23.2 

44.6 

47.5 

33.9 

37.2 

EKZ/188/C 1 N 

20.  1 

21.3 

13.6 

15.5 

27.2 

29.  3 

15.9 

18.6 

EKZ/041/05N 

25.  3 

25.  5 

24.3 

24.  7 

54.9 

57.7 

42.5 

45.6 

EKZ/364/C6N 

25.  5 

25.6 

24.7 

24 . 9 

58.5 

60.3 

47.  1 

48.9 

EKZ/349/07N 

22.0 

22.8 

16.  5 

18.7 

31.6 

34.  2 

20.4 

23.8 

EKZ/333/06N 

25.  1 

25.4 

23.4 

24.  C 

48.6 

51.2 

37.  1 

4C.0 

EKZ/294/06N 

24.8 

25. C 

22.4 

23.  3 

46.6 

48.5 

34.  1 

37.  1 

EKZ/159/C  IN 

24.8 

25.  2 

23. C 

23.  9 

45.3 

48.7 

36.0 

39.  1 

NFV/230/14N 

21.8 

22.6 

17. C 

19.4 

32. 1 

34.9 

21.8 

26.4 

CAN/310/22N 

25.  4 

25.  4 

24.3 

24.4 

55.0 

56.4 

44.  1 

45.2 

NZH/270/06N 

25.5 

25.5 

24.9 

25.0 

62.7 

63.2 

50.9 

51.6 

TAD/189/04N 

19.3 

19.9 

12.3 

13.4 

24.5 

25.6 

13.9 

15.5 

IRA/102/19N1 

24.7 

24.8 

23.2 

23.5 

4Q.9 

50.8 

41.2 

42.3 

KUR/077/03N 

20.  1 

20.3 

14.  4 

15.7 

29.5 

30.9 

18.0 

20.2 

KAZ/181/04N 

23.9 

24.3 

21.3 

22.3 

43.  9 

47.4 

33.2 

36.  2 

KH1/073/12N 

25.  1 

25.2 

23.7 

24.2 

51.1 

53.4 

38.1 

4''. 7 

KAZ/356/06N 

25.  6 

25.6 

24.9 

24.9 

56.8 

56.  9 

47.  1 

47.4 

KAZ/282/06N 

25.0 

25.1 

23.0 

23.  6 

46.2 

48.4 

34.7 

37.  1 

0RA/191/16N 

25.  2 

25.3 

23.9 

24.  2 

50.5 

51.8 

40.  9 

42.4 

KAZ/157/04N 

25.3 

25.4 

24.2 

24.5 

50.7 

52.9 

40.4 

42.7 

KAZ/145/04N 

23.  9 

24.4 

20.9 

21.9 

39.4 

4 1.2 

28.5 

31. C 

KAZ/1 15/03N 

25.  6 

25.6 

24.8 

24.9 

59.2 

59.6 

47.0 

48.0 

KAZ/170/04N 

25.4 

25  4 

24.0 

24.  3 

52.3 

53.  6 

39.6 

41.6 

KAZ/081/04N 

25.5 

25.5 

24.4 

24.7 

54.1 

55. 2 

42.7 

44.6 

WRS/295/05N 

24.7 

24.8 

22.  7 

22.9 

51.2 

51.  1 

42.0 

43  . 1 

TUG/053/21 N 

14.  8 

14.3 

7.2 

7.  3 

21.3 

21.2 

9.6 

10. C 

ITA/039/12N 

12.6 

12.0 

4.8 

4 . 8 

20.2 

20.4 

8.9 

9.6 

ITA/036/C7N 

13.7 

14.3 

8.  1 

9.  1 

20.  9 

22.  3 

9.8 

11.4 

ITA/C36/05N 

17.6 

18. C 

9.2 

1 . 1 

23.7 

24.  3 

1C.  6 

11.7 

ITA/036/03N 

1 3.  2 

13.5 

3.2 

3.5 

16.4 

16.4 

5.0 

5.3 

ITA/036/0  IN 

20.  2 

20.8 

14.  2 

15.  1 

31.9 

32.  1 

18.3 

19.7 

GPE/074/15N 

24.6 

24.9 

21.6 

22.5 

43.6 

44.8 

3C  .6 

33.4 

YUG/100/02N 

19.6 

20.7 

12.  3 

14.  1 

27.4 

30.  1 

18.2 

19.7 
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TABLE  III-2 


DETECTOR  OUTPUT  (dB)  (1.  5-2.  5 Hz  BAND) 
(PAGE  2 OF  4) 


Coh.  Incoh. 


Coh. 


Incoh. 


Event  I.  D. 

A.  D. 

D.S. 

A.E 

TUR/067/23N 

13.  5 

13.8 

7.  1 

GRE/109/02N 

21.0 

21.7 

17.  6 

AUS/169/09N 

20.  9 

21.2 

12.  9 

ITA/166/18N 

17.2 

18.0 

12.7 

GRE/167/00N 

23.0 

2 3.2 

19.2 

ITA/173/15N 

10.9 

1C. 4 

5.9 

YUG/176/07N 

19.7 

20.  1 

13. C 

YUG/177/04N 

21.3 

21.3 

13.8 

YUG/180/01N 

13.4 

12.8 

6.4 

GRE/187/18N 

18.6 

18.5 

12.  6 

WES/262/11N 

22.6 

22.9 

19.  1 

TDR/161/09N 

20.9 

21.2 

18.  1 

TUR/126/04N 

17.4 

18.0 

13.8 

BLS/210/19N 

19.5 

19.8 

13. C 

TUR/276/07N 

16.3 

17.0 

10.5 

TUR/278/18N 

18.2 

18.8 

10.5 

GRE/207/P1N 

22.6 

23. C 

18.  3 

AEG/212/01N 

16.4 

16.8 

12.7 

MRS/191/07N 

22.4 

22.3 

18.  1 

MRS/248/C7N 

23.  2 

23.8 

19.7 

ITA/035/32N 

16.  6 

17.2 

11.2 

ITA/035/04N 

12.7 

13.0 

2.4 

ITA/035/09N 

13.  2 

13.3 

9.2 

ITA/035/19N 

13.0 

13.0 

4.5 

ITA/035/17N 

19.0 

19.6 

11.9 

ADR/037/01N 

21.9 

22.4 

16.  2 

GRE/200/13N 

19.2 

19.5 

13.3 

ITA/018/23N 

17.  6 

18.4 

10.5 

YDG/067/05N 

17.8 

18.2 

10. C 

YUG/052/23N 

18.  7 

19.2 

11.2 

CAU/079/03N 

12.8 

13.  1 

4.4 

SWR/076/02N 

13.  1 

12.6 

2.  7 

HRS/202/15N 

11.5 

10.9 

1.7 

SWR/205/1 IN 

13.9 

14.3 

7.5 

TUR/173/05N 

11.5 

10.7 

5.4 

TUR/167/14N 

16.  1 

15.  6 

8.  1 

GRE/175/07N 

13.  6 

14.5 

6.4 

TUR/175/04N 

14.0 

14.7 

4.9 

WRS/179/12N 

14.  6 

14.5 

6.6 

RDM/251/04N 

19.0 

19.5 

11.  3 

GRE/153/13N 

17.  8 

17.9 

9.C 

GRE/157/10N 

22.  4 

22.7 

17.  9 

CRI/204/C5N 

24.9 

25.  1 

23.4 

TUP/212/19N 

16.  1 

15.8 

7.8 

AIFS/0GS/04n 

15.  3 

14.9 

6.  3 

RUL/068/22N 

16.  1 

17.5 

7.4 

1.  D.  S. 

A.D. 

D.S. 

A.D 

. D.S 

7.  8 

16.3 

17.4 

7.8 

8.7 

18.0 

36.7 

37.9 

25.  1 

2 6.  5 

13.8 

29.5 

30.  1 

15.7 

16.9 

13.  1 

27.0 

28.9 

20.3 

22.  1 

19.  7 

45.  1 

46.5 

33.8 

35.  1 

5.9 

12.3 

11.5 

6 . 6 

6. 

13.7 

26.6 

27.2 

15.7 

16.7 

15.  0 

28.0 

28.4 

16.2 

17.9 

6.6 

14.9 

14.2 

9.  1 

8.9 

13.  3 

27.7 

28.'^ 

16.4 

17.7 

20.  1 

40.0 

40.6 

33.5 

34.9 

18.  8 

32.  3 

34.0 

26.2 

27.5 

14.3 

29.  1 

30.4 

21.8 

23.  1 

14,3 

25.  3 

25.8 

15.0 

16.7 

11.1 

21.7 

22.8 

12.4 

13.1 

11.5 

23.7 

24.0 

11.9 

13.2 

18.9 

36.4 

37.6 

25.  9 

27.3 

13.3 

32.  9 

33.6 

2C.7 

22.  ^ 

19.  1 

37.5 

40.  1 

37.  1 

38.6 

21.0 

45.9 

47.1 

34.  3 

36.  7 

11.7 

28.7 

29.3 

16.9 

18.2 

2.  9 

15.6 

15.  8 

2.7 

3.2 

9.7 

23.5 

23.7 

12.7 

13.6 

4.9 

19.7 

20.0 

7.8 

8.  3 

12.  9 

26.6 

27.  1 

14.6 

16.0 

17.  1 

34 . 4 

35. C 

22.5 

23.8 

14.  0 

29.0 

30.3 

17.3 

18.7 

12.  5 

25.  3' 

26.8 

12.7 

15.4 

10.8 

24.0 

24.  8 

11.6 

12.6 

12.  4 

26.7 

2®.0 

13.6 

15.2 

4.  3 

16.4 

17.  2 

5.9 

6.  1 

2.9 

15.0 

14.7 

2.7 

3.0 

2.C 

13.2 

12.9 

2.0 

2.3 

8.0 

2C.4 

20.6 

10.6 

11.3 

5.  8 

15.6 

14.8 

7.5 

7.9 

8.7 

25.2 

25.  8 

13.6 

14.9 

8.  2 

19.9 

21.6 

9.  1 

11.8 

6.2 

16.6 

17.6 

5.2 

6.  5 

7.0 

17.4 

17.0 

7.^ 

8.2 

12.  C 

26.0 

26.3 

14.9 

15.8 

9.9 

20.5 

2'^. 9 

10.0 

1C. 9 

18.  6 

33.  7 

35.4 

22. « 

24. 2 

23.7 

51.2 

51.5 

W'  .2 

41.1 

8.  1 

19.7 

19. 2 

9.4 

8.6 

6.  8 

19.5 

19.6 

7.5 

8. : 

9.3 

2(  . 2 

21.7 

8.  1 

1C.  4 
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TABLE  III-2 


DETECTOR  OUTPUT  (dB)  (3.  0-4.  0 Hz  BAND) 

(PAGE  3 OF  4) 

Coh.  Incoh.  Coh.  Incoh. 


Event  I.  D. 

A.D. 

, D.S. 

A.D, 

KUR/1U2/12N 

8.0 

5.8 

1.2 

KAZ/070/0UN 

24.0 

23.8 

22.2 

EK2/3U5/0UN 

24.6 

24.9 

23.6 

EK1/3U5/0UN 

24.1 

24.5 

23.2 

SWR/277/08N 

24.2 

24.0 

21.4 

NZM/2U 1/05N 

25.  4 

25.  4 

24.  1 

EKZ/363/0UN 

20.7 

20.6 

15.0 

EKZ/229/03N 

2 3.  6 

23.8 

21.5 

EKZ/1B8/01N 

18.  1 

18.5 

11.7 

EKZ/04 1/05N 

24.9 

25.0 

23.3 

EKZ/2U6/C8N 

23.4 

23.7 

20.  1 

EKZ/239/03N 

24.7 

2 5.1 

23.5 

EKZ/36U/06N 

24.8 

24.8 

23.  3 

EKZ/349/07N 

21.8 

22.3 

16.  1 

EKZ/333/06N 

24.4 

24.6 

22.6 

EKZ/294/06N 

24.2 

24.4 

22.  1 

EKZ/159/0 IN 

24.4 

24.6 

21.8 

WRS/277/10N 

22.4 

22.  6 

21.7 

NEV/230/14N 

5.2 

7.4 

1.9 

KAZ/181/04N 

24.5 

25.0 

21.7 

nRA/C82/06N 

24.3 

24.  1 

22.8 

CAN/310/22N 

25.  4 

25.  4 

24.9 

NZM/270/06N 

25.2 

25.2 

24.7 

TAD/188/04N 

19.0 

20.3 

13.7 

IRA/102/19N1 

22.6 

23.2 

20.8 

KUR/077/03N 

10.  4 

12.7 

8.  1 

KM1/073/12N 

23.  P 

24.3 

21.9 

JAP/063/07N 

4.7 

4.3 

3.6 

KAZ/356/06N 

24.4 

24.7 

24.3 

KAZ/145/04N 

23.8 

24.2 

22.4 

KAZ/1 15/03N 

24.  6 

24.7 

24.8 

KAZ/170/04N 

24.3 

24.5 

23.2 

KAZ/08 1/04N 

24.8 

24.9 

24.5 

KAZ/282/06N 

24.4 

24.8 

23.3 

URA/191/lf-N 

24.  1 

24.3 

22.  6 

KAZ/157/04N 

23.7 

24  . U 

23.5 

WRS/245/05N 

24.3 

24.4 

2 3.3 

YUG/063/21N 

13.0 

13.2 

4.  1 

ITA/036/07N 

13.0 

12.7 

6.  1 

TTA/036/05N 

9.7 

9.3 

2.  6 

ITA/036/03N 

8.  1 

7.8 

1.5 

TTA/036/0 IN 

17.8 

18.2 

10. 9 

TTA/039/12N 

9.6 

4.3 

2.8 

ADR/037/0  IN 

20.8 

21.4 

1 4.  4 

WES/262/1  IN 

23.4 

23.  1 

21.6 

D.S. 

A.  D. 

D.S, 

A.  D 

. D.S 

2.3 

8.4 

6.  2 

1.3 

2.  1 

22.5 

4 9.6 

50.2 

38.9 

40.  1 

24.1 

59.7 

61.2 

48.6 

50.6 

23. 7 

60.4 

62.2 

48.5 

50.3 

21.5 

47.0 

47.7 

38.9 

39.9 

24.  2 

57.3 

57.7 

52.7 

52.9 

15.9 

27.5 

25.8 

18.2 

18.  1 

22.2 

42.8 

44.  1 

3 3.0 

34.8 

12.7 

23.4 

24.0 

14.3 

15.6 

23.  9 

51.3 

54.5 

40.1 

43.5 

20.  8 

38.6 

39.5 

27.9 

29.8 

24.  2 

51.8 

54.8 

41.9 

44.5 

23.6 

51.6 

52.8 

4 1.2 

43.3 

17.  5 

33.  1 

34.7 

20.8 

23.0 

23.  1 

47.5 

48.9 

36.2 

38.2 

22.6 

45.8 

45.7 

34.6 

36.  1 

22.6 

40.5 

44.3 

31.7 

34.5 

22.  1 

3«.  8 

4 1.0 

38.0 

38.5 

3.  1 

6.  1 

8.  1 

2.  1 

3.4 

23. C 

44.2 

49.4 

33.  9 

38.0 

22.8 

53.3 

53.8 

45.9 

46.5 

25.0 

55.5 

55.  2 

47.7 

48.7 

24.8 

58.1 

59.0 

51 . 4 

52.2 

15.  4 

24.6 

27.  1 

16.2 

18.6 

21.2 

3 8.1 

41.5 

40.2 

42.2 

9.5 

11.8 

14.8 

9.5 

11.1 

23.  1 

41.7 

46.  9 

32.8 

38.  3 

3.9 

7.7 

9.7 

5.  1 

5.9 

24.4 

47.0 

52.  2 

48.6 

50. 1 

22.8 

40.1 

42.4 

33.9 

35.5 

24.9 

55.  2 

56.8 

51.4 

52.8 

23.6 

45.5 

46.7 

39.0 

40.5 

24.7 

52.5 

56.3 

47.6 

49.6 

23.8 

46.2 

49.  R 

38.4 

41.0 

22.9 

42.8 

44.7 

38.5 

39.  1 

24.  1 

42.2 

46.  3 

39.  1 

42.  0 

23.3 

47.5 

48.6 

4 1.6 

42.9 

4.4 

15.7 

16.2 

5.2 

5.6 

6.4 

17.0 

16.8 

7.5 

7.8 

2.  9 

12.8 

11.7 

3.6 

4.0 

1.9 

8.6 

8.  3 

2.2 

2.7 

11.9 

2 3.5 

23.5 

12.4 

13.7 

2.  9 

12.2 

11.8 

3. -7 

3.8 

16. C 

28.8 

31.1 

17.7 

20.  3 

21.5 

4 8. 5 

48.6 

43.8 

4 3.5 
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table  I1I-2 

DETECTOR  OUTPUT  (dB)  (3.  0-4.  0 Hz  BAND) 
(PAGE  4 OF  4) 


Event  I.  D. 
TUR/161/09N 
TUR/126/04N 
BLS/210/19N 
ITA/035/19N 
ITA/035/17N 
ITA/035/09N 
ITA/035/0UN 
ITA/035/02N 
GRE/07U/15N 
YnG/100/02N 
T0R/067/23N 
GRE/109/02N 
ITA/166/18N 
GRE/167/00N 
YTIG/180/01N 
GRE/187/18N 
YUG/177/04N 
YOG/176/07N 
AUS/169/09N 
ITA/173/15N 
TUR/278/18N 
TUR/276/07N 
GRE/207/01N 
WRS/191/07N 
HRS/248/07N 
AEG/212/01N 
GRE/200/13N 
ITA/018/23N 
YUG/067/05N 
YOG/052/23N 
BOL/068/22N 
CAU/079/03N 
Rtin/251/04N 
A0S/005/04N 
SWR/076/02N 
HRS/202/15N 
SHR/205/1 IN 
TUR/173/05N 
T0R/167/14N 
WRS/179/12N 
TUR/175/04N 
GRE/175/07N 
GRE/153/13N 
GRE/157/10N 
CRI/204/05N 
TUR/21 2/19N 


A.  D^'^rs. 

A.  ifP’TSh. 

4.9  7.7 

3.9 

4.4 

5.3  7.1 

9.3 

9.4 

4.6  6.4 

4.0 

4.2 

16.2  15.5 

1 

7.9 

15.6  14.7 

6.2 

6.8 

11.6  11.4 

4.6 

4.  1 

7.8  7.3 

1.5 

1.5 

15.9  15.9 

11.8 

12.5 

18.7  20.9 
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SECTION  IV 

RESULTS  AND  DISCUSSIONS 


A.  INTRODUCTION 

This  section  presents  the  results  of  the  computation  and  proces- 
sing of  the  events  and  noise  samples.  Subsection  B illustrates  the  computation 
of  beams,  and  the  associated  averages  (or  envelopes)  and  the  detector  output. 
Subsection  C discusses  the  false  alarm  probability  measurements.  Following 
in  Subsection  D is  the  discussion  of  detection  probability  measurements.  The 
operating  characteristics  of  the  detectors  are  presented  and  compared  in  the 
last  subsection. 


B.  ILLUSTRATION  OF  COMPUTATION 

Figure  IV-1  shows  coherent  and  incoherent  adjusted-delay  and 

diversity-stack  beams  (AB),  and  the  associated  short-term  averages  (STA)  and 

long-term  averages  (LTA)  for  a signal  of  a presumed  explosion  from  Kazakh. 

In  this  figure,  all  array  beams  a_*e  plotted  with  the  same  scale  factor  but  the 

associated  averages  are  enlarged  by  a factor  l.Z  x (AB)  STA)  in  each 

max  max 

case.  The  shape  of  the  STA  approximately  reflects  the  envelope  and  duration  of  the 
P-wave  although  it  is  delayed  by  the  1.5  second  integration  time.  An  adjusted- 
delay  beamformer  is  used  in  the  NORSAR  on-line  system,  but  the  diversity- 
stack  beamformer  yields  the  better  signal-to-noise  processing  gain,  and  hence, 
potentially  provides  an  improvement  in  detectability.  Figure  IV-2  illustrates 
the  same  computations  for  a noise  sample  from  day  214  of  1971.  Because  the 
associated  STA  and  LTA  averages  of  the  noise  are  also  multiplied  by  the  factor 
given  above,  STA  and  LTA  appear  slightly  larger  than  AB. 
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ILLUSTRATION  OF  ARRAY  BEAMS,  STAS,  LTAS 
(DAY  214  OF  1971  NOISE,  1.5-2.  5 Hz) 


Figure  IV- 3 shows  the  beams  ana  the  associated  envelopes  and 
long-term  averages  for  the  same  event  shown  in  Figure  IV- 1,  but  where  ;he 
Hilbert  transform  is  used  for  the  envelope  formation.  All  the  traces,  beams 
and  envelopes,  have  the  same  plot  scale  factor  in  this  figure.  This  method 
produces  an  obviously  more  accurate  representation  of  the  envelope  than  the 
STA  technique  and  greatly  increases  the  detectability  of  P-wave  amplitudes  as 
shown  in  the  figure.  Figure  IV-4  repeats  the  same  illustration  for  a noise 
sample  from  day  214  of  1971.  Note  that  due  to  the  absence  of  noise  cancella- 
tion involved  in  incoherent  array  beamforming,  the  incoherent  envelopes  have 
the  higher  amplitudes;  but  that  the  variability  of  the  STA  output  is  less  for  the 
incoherent  beam. 

Figures  IV-5  and  IV-6  illustrate  the  detector  output  for  an  event 
from  Kazakh  and  a noise  sample  from  day  214  of  1971,  respectively.  In  these 
figures,  the  output  is  the  running  signal-to-noise  ratio,  SNR(t)  of  equation 
(11-13).  Since  each  trace  has  a different  plot  scale-factor,  it  does  not  reflect 
the  performance  of  detectors.  For  Figure  IV-5,  the  peak  values,  (STA/ LTA) 

max 

are  shown  in  Table  III-2.  Figures  IV-7  and  IV-8  show  the  same  illustrations 
for  the  event  in  Figure  IV-5  and  the  noise  in  Figure  IV-6,  but  using  a Hilbert 
transformation.  The  higher  variabilities  of  detector  output  by  use  of  Hilbert 
transformation  reflect  more  degrees  of  freedom  for  the  analytic  envelopes  than 
for  the  STA  envelopes. 

C.  THE  MEASURED  FALSE  ALARM  PROBABILITY 

As  an  illustration.  Figure  IV-9  shows  a histogram  of  detector 
output  from  an  incoherent  noise  beam.  This  curve  was  normalized  to  one  at 
maximum  density.  The  measured  mean  value  and  standard  deviation  are  -0.  053 
dB  and  0.  691  dB,  respectively.  In  the  lower  part  of  the  figure  is  shown  the 
number  of  counts  within  a certain  output  range.  For  example,  there  are  9 
counts  (out  of  a total  of  2900)  in  the  range  from  2.  0 to  2.  2 dB. 
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FIGURE  IV-3 

AN  ILLUSTRATION  OF  ARRAY  BEAMS, 
ENVELOPES  AND  LTAS  USING  HILBERT  TRANSFORM 
(KAZ/282/06N,  1.5-2.  5 Hz) 


AN  ILLUSTRATION  OF  ARRAY  BEAMS,  ENVELOPES,  AND  LTAS  USING  HILBERT  TRANSFORM 

(DAY  214  OF  1971  NOISE,  1.  5-2.  5 Hz) 


Coherent  Diversity-Stack  Beamforming 


FIGURE  IV-5 

AN  ILLUSTRATION  OF  DETECTOR  OUTPUT,  SNR 
(KAZ/282/06N,  1.5-2. 5 Hz) 
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Coherent  Adjusted-Delay  Beamforming 


Incoherent  Adjusted-Delay  Beamforming 


Coherent  Diversity-Stack  Beamforming 


Incoherent  Diversity-Stack  Beamforming 


FIGURE  IV-7 


AN  ILLUSTRATION  OF  DETECTOR  OUTPUT  USING  HILBERT  TRANSFORM 

(KAZ/282/06N,  1.5-2. 5 Hz) 
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Figures  IV-10  and  IV-11  show  the  false  alarm  probability  dis- 
tributions for  the  coherent  and  incoherent  detectors  computed  from  ten  noise 
samples  for  the  adjusted-delay  beams  in  the  1.5-2.  5 Hz  and  3.  0-4.  0 Hz  pass- 
bands,  respectively.  The  standard  deviations  of  the  probability  for  a given 
detector  output  level  among  the  ten  noise  samples  are  indicated  for  several 
discrete  points.  The  arrows  indicate  that  the  standard  deviation  is  greater 
than  the  mean  value  for  that  specific  level,  which  means  that  the  experimental 
false  alarm  probabilities  can  possibly  be  zero.  Figures  lV-12  and  lV-13  show 
the  false  alarm  probability  distributions  for  the  diversity-stack  beams  in  the 
1»5-2.5  Hz  and  3.  0-4.  0 Hz  passbands,  respectively.  They  were  computed 
from  the  same  ten  noise  samples. 

The  noise  beams  were  formed  at  infinite  velocity  (vertical  beams) 
since,  as  will  be  discussed  later,  the  noise  is  essentially  uncorrelated  between 
subarrays  at  NORSAR  (Barnard  and  Whitelaw,  1972),  For  the  diversity-stack 
beams,  the  probability  curves  closely  follow  a Gaussian  distribution  but,  the 
tails  of  the  adjusted-delay  curves  behaved  erratically  and  deviated  from  a Gaus- 
sian distribution  in  both  passbands.  This  is  partially  due  to  the  limited  data 
base.  The  diversity-stack  curves  are  more  well-behaved  because  the  weighting 
coefficient  of  each  channel  was  proportional  to  the  reciprocal  of  the  noise  pow- 
er and  thus  tended  to  suppress  the  erratic  channels. 

The  discussion  of  false  alarm  rate  measurement  will  be  based  on 
the  diversity-stack  beams  because  it  is  more  uniformly  distributed  and  it  is  be- 
lieved that  the  results  will  be  more  reliable.  Comparing  the  results  between 
Figures  IV-12  and  lV-13,  the  false  alarm  probability  in  the  1.5-2.  5 Hz  pass- 
band  is  higher  than  that  in  the  3,  0-4.  0 Hz  passband  for  the  coherent  detector 
while  the  opposite  is  true  for  the  incoherent  detector.  For  example,  having  the 

_3 

false  alarm  probability  at  the  10  level,  the  decision  threshold  is  7.5  dB  in  the 
1,  5-2.  5 Hz  passband  and  7.0  dB  in  the  3.  0-4.  0 Hz  passband  for  the  coherent  de- 
tector. This  difference  increases  at  lower  probability  levels.  For  an  incoherent 
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FALSE  ALARM  PROBABILITY  FOR  ADJUSTED-DELAY  BEAM 
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detector  having  the  same  10  level,  the  decision  threshold  is  2,3  dB  in  the 
1.  5-2.  5 Hz  pas'jband  and  3.  5 dB  in  the  3.  0-4.  0 Hz  passband.  Therefore,  the 
false  alarm  rate  is  much  higher  in  the  higher  frequency  passband  for  an  in- 
coherent detector. 

The  explanation  of  this  result  is  given  below.  Previous  short- 
period  analysis  indicated  that  the  noise  was  almost  completely  uncorrelated 
among  subarrays  because  of  the  large  distances  involved  (IBM,  1970).  The 
multiple  coherence  between  the  reference  sensor  of  subarray  1 and  the  refer- 
ence sensors  of  subarrays  2,  3,  4,  6,  and  8 for  one  noise  sample  was  about 
0.  14  in  the  1. 5-2.  5 Hz  passband  and  0.  1 in  the  3.  0-4.  0 Hz  passband  (Barnard 
and  Whitelaw,  1972,  Figure  11-9).  This  could  account  for  the  higher  false  alarm 
probability  in  the  1.5-2.  5 Hz  passband  than  in  the  3,  0-4.  0 Hz  passband  for  the 
coherent  detector.  However,  the  noise  coherence  was  higher  in  the  3.  0-4.  0 Hz 
passband  than  in  the  1,5-2.  5 Hz  passband  within  a subarray  at  NORSAR  (Bar- 
nard and  Whitelaw,  1972,  Figure  11-8;  Capon,  1972,  Figure  IV-15).  This  could 
account  for  the  greater  false  alarm  rate  in  the  3.  0-4.  0 Hz  passband  than  in  the 
1.5-2.  5 Hz  passband  with  the  incoherent  beam  because  it  is  essentially  the  aver- 
age of  the  subarray  beams. 

Lacoss  (1972)  using  NORSAR  data  theoretically  estimated  the 

-5 

false  alarm  probability  for  the  10.  5 dB  level  at  1.1  x 10  for  the  coherent  de- 
tector. From  Figures  IV-12  and  IV-13,  the  decision  threshold  would  be  slight- 
ly lower  than  that  for  the  coherent  detector.  But,  for  the  incoherent  detector, 
the  threshold  must  be  higher  than  the  3.  0 dB  level  in  the  1.  5-2.  5 Hz  passband 
and  the  4.5  dB  level  in  the  3,  0-4.0  Hz  passband.  Table  IV - 1 shows  the  statis- 
tical parameters  derived  from  the  observed  probability  curves  shown  in  Figures 
IV- 10  through  IV-13.  In  using  these  parameters  to  estimate  the  false  alarm 
probability  it  must  be  noted  that  the  estimated  probability  could  be  higher  than 
actually  shown  in  Figures  IV-10  through  IV-13  in  the  tail  because  the  experi- 
mental curves  go  to  zero  there.  For  example,  with  the  diversity-stack  coherent 
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detector,  in  order  to  achieve  a false  alarm  rate  less  than  3.  17  x 10  , we 

must  have  (S  + 0.  38)  / 2, 70  = 4.  0 or  S = 1 0.  42  dB.  This  is  higher  than  what 


Figure  IV-13  indicates. 


D.  DETECTION  PERFORMANCE 


Comparison  of  Outputs  From  Coherent  and  Incoherent  Detectors 


One  of  the  objectives  of  this  study  was  to  compare  the  perfor- 
mance of  coherent  and  incoherent  detectors  on  signals  from  close  events  where 
conventional  coherent  beamforming  is  occasionally  ineffective.  This  problem 
is  compounded  by  signals  with  particularly  strong  high-frequency  energy  such 
as  from  close  events  where  attenuation  of  that  energy  is  not  as  great.  The  an- 


alysis of  the  91  events  listed  in  Table  III-l  was  oriented  toward  a study  of  these 


effects. 


Figures  IV-14  through  IV -21  are  the  plots  of  the  coherent  detec- 
tor output  (in  dB)  versus  the  incoherent  detector  output.  Tin-  first  four  figures 
illustrate  the  results  in  the  1,5-2. 5 Hz  passband  using  adjusted-delay  and 


diversity-stack  beams  for  the  (STA/LTA)  mode  and  the  (ST A)  /LTA 

max  max 


mode,  respectively.  The  next  four  figures  repeat  the  same  illustrations  for 
the  3.  0-4.  0 Hz  passband.  The  presumed  explosions  and  the  earthquakes  are 
identified  by  different  symbols  as  indicated  on  the  figures. 


The  coherent  detector  performed  better  than  the  incoherent  de- 
tector in  the  lower  frequency  band  than  in  the  higher  frequency  band.  Figures 
for  the  3.  0-4.0  Hz  passband  seem  to  suggest  that  the  incoherent  detector  tends 


to  perform  better,  particularly  for  the  lower  magnitude  events.  Hence,  for  de- 
tecting the  small  events,  the  incoherent  detector  operating  in  the  higher  frequency 
band  would  be  valuable.  However,  for  better  comparison  of  the  detector's  per- 


formance, it  is  necessary  to  compare  the  operating  characteristics  of  the  de- 


tectors involved. 
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FIGURE  IV- 14 

COHERENT  ADJUSTED-DELAY  DETECTOR  OUTPUT  VERSUS 
INCOHERENT  DETECTOR  OUTPUT 
((STA/LTA)^^  IN  THE  1.5-2.  5 Hz  PASSBAND) 
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Coherent  Detector  Output  (dB) 

FIGURE  IV-15 

COHERENT  DIVERSITY-STACK  DETECTOR  OUTPUT  VERSUS 

INCOHERENT  DETECTOR  OUTPUT 

( (STA/LTA)  IN  THE  1.  5-2.  5 Hz  PASSBAND) 
max 
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FIGURE  IV-16 


COHERENT  ADJUSTED-DELAY  DETECTOR  OUTPUT  VERSUS 
INCOHERENT  DETECTOR  OUTPUT 
( (STA)  /LTA  IN  THE  1.  5-2.  5 Hz  PASSBAND) 
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To  summarize  the  detector  outputs,  Table  IV-2  shows  the  aver- 
age output  for  the  coherent  and  incoherent  detectors.  In  the  updated  mode, 
i,  e.  , (STA/LTA)  , the  coherent  detector  was  2-3  dB  better  than  the  in- 
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coherent  detector  in  the  two  passbands  for  the  presumed  explosions.  For  the 
earthquake  ensemble  it  was  5-6  dB  better.  The  presumed  explosion  group 
had  higher  signal-to-noise  ratios  and  resulted  in  smaller  signal  loss  in  array 
beamforming  because  picking  the  correct  delay  anomalies  was  easier  for  high- 
er SNR  signals.  For  example,  for  the  adjusted  delay  detector,  the  average 
decrease  from  the  presumed  explosion  group  to  the  near-regional  earthquakes 
was,  for  the  coherent  detector,  6.2  dB  in  the  1.5-2.  5 Hz  passband  and  8.3  dB 
in  the  3.  0-4.  0 Hz  passband,  respectively,  and  9.  9 dB  in  the  1. 5-2.  5 Hz  pass- 
band  and  12.0  dB  in  the  3.  0-4.0  Hz  passband,  respectively,  for  the  incoherent 
detector. 


For  the  frozen  and  averaged  LTA  mode,  using  (STA)  /LTA, 
the  coherent  detector  was  10  - 11  dB  better  than  the  incoherent  detector  in  the 
1.5-2.  5 Hz  passband  and  7 - 8 dB  better  in  the  3.  0-4.  0 Hz  passband.  In  this 
mode,  the  SNR  estimate  is  more  physically  meaningful.  The  theoretical  array 
gain  for  22  channels  is  13.4  dB  and  the  signal  loss  in  beamforming  at  NORSAR 
was  about  3 - 4 dB  in  the  signal  band  (Ringdal  and  Whitelaw,  1973).  The  small- 
er improvement  of  the  coherent  detector  over  the  incoherent  detector  in  the 
3.0-4.  0 Hz  passband  reflects  the  higher  signal  loss  due  to  beamforming  in  the 
higher  frequency  passband. 


Envelope  detectors  using  Hilbert  transformation  were  applied  to 
four  events  - two  presumed  explosions  and  two  earthquakes.  Because  each 
data  point  E(t)  in  either  the  coherent  or  the  incoherent  beamforming  process 
represents  an  independent  estimate  of  amplitude,  the  theoretical  maximum  in- 
crease from  STA-envelope  detectors  for  coherent  beamforming  design  would  be 
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coherent  detector  and  from  0.  5 dB  to  8.6  dB  for  the  incoherent  detector.  The 
average  increase  was  7.  3 dB  and  4.  1 dB  for  the  (STA/ LTA)  and 
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(STA)  / LTA  modes,  respectively,  with  the  coherent  diversity-stack  detec- 
max  

tor  and  5.  1 dB  and  2.  1 dB  for  the  (STA/ LTA)  and  (STA)  /LTA  modes 
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with  the  incoherent  diversity-stack  detector.  Thus  the  operational  mode  de- 
tector registered  a more  significant  increase  than  the  frozen  and  averaged 
LTA  mode  detector. 

2.  A Simulation  Study 

In  order  to  determine  the  detector  performance  at  different 
event  magnitudes,  a simulation  study  was  conducted  by  burying  a signal  scaled 
over  a range  of  amplitudes  in  subarray  beams  of  noise  and  subsequently  gen- 
erating the  coherent  and  incoherent  beams  for  the  detection  system.  The  sub- 
array beam  amplitudes  were  obtained  by  computing  the  ratio  of  the  signal  peak 
amplitude  to  the  noise  RMS  value  for  each  subarray  and  averaging  the  ratios 
over  the  full  array.  These  ratios,  which  are  expressed  in  decibels,  are  call- 
ed the  subarray  SNR's. 

Figure  IV-22  shows  the  detector  output  versus  the  mean  sub- 
array signal-to-noise  ratio  for  the  presumed  explosion  KAZ/282/06N.  Curves 
Aj  and  A^  are  for  the  frozen  and  averaged  long-term  average  mode,  i.  e.  , 

(STA)  /LTA  , of  coherent  detector  operating  in  the  1.  5-2.5  Hz  and  3.0-4.  0 
max  

Hz  passbands  , respectively.  Curves  B and  B are  for  the  (STA)  / LTA 

mode  of  incoherent  detector  output  in  the  1. 5-2.  5 Hz  and  3.  0-4.  0 Hz  pass- 

bands,  respectively.  Curves  C^  and  C^  are  for  the  updated  long-term  aver- 

ave  mode,  i.  e.  , (STA/ LTA)  , of  the  coherent  detector  output  in  the  1.5- 

max 

2.  5 Hz  and  3.0-4,  0 Hz  passbands,  respectively,  while  curves  and  are 
for  the  updated  LTA  mode  of  incoherent  detector  in  the  1.  5-2.  5 Hz  and  3.  0- 
4.  0 Hz  passbands,  respectively. 
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FIGURE  IV-22 

MEAN  SUBARRAY  SNR  VERSUS  DETECTOR  OUTPUT 

(KAZ/282/06N) 


For  a given  mean  subarray  SNR,  the  incoherent  detector  per- 
formed better  in  the  3.  0-4.  0 Hz  passband  while  the  reverse  was  true  foi-  the 
co^'erent  detector.  This  reflects  the  fact  that  the  signal  loss  was  higher  in 
the  3.  0-4.  0 Hz  passband  than  in  the  1.  5-2.  5 Hz  passband  for  the  coherent  de- 
tector. Typical  SNR  detection  thresholds  used  at  NORSAR  are  11  dB  for  the 
coherent  detector  and  4 dB  for  the  incoherent  detector  (Ringdal  et  al.  , 1972). 
Using  these  levels  as  the  basis  for  defining  the  false  alarm  rates  in  order  to 
compare  the  detector's  performance,  then,  as  the  mean  SNR  decreases,  the 
detector  output  reaches  these  levels  in  the  order:  curves  C^.  D . B . A . D . 
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B^,  Cj,  and  A^.  Direct  interpretation  of  this  order  suggests  that  curve  A^, 
the  (STA)  /LTA  mode  of  coh  erent  detector  operating  in  the  1.  5-2.  5 Hz 

ITlaX 

passband,  yields  best  performance  and  the  (STA/ LTA)  mode  of  coherent 

max 

detector  operating  in  the  3.  0-4.  0 Hz  passband  ranks  last  in  performance. 


For  the  overall  performance  of  the  two  groups  of  curves,  the 
incoherent  detector  has  better  performance  than  the  coherent  detector  in  the 
3.  0-4.  0 Hz  passband,  but  not  in  the  1.  5-2.  5 Hz  passband.  Also,  for  a given 
mean  SNR,  the  difference  in  detector  output  between  the  two  passbands  is 
about  2 dB  for  the  coherent  detector,  and  is  less  than  1 dB  for  the  incoherent 
detector. 


Figure  lV-23  shows  the  results  of  the  same  simulation  study 
for  the  event  WRS/248/07N.  which  was  11.9  degrees  away  and  7 km  deep. 
Curves  A^,  A^,  B^,  B^,  C^,  C^,  D^,  and  correspond  to  the  same  cases 
as  in  Figure  IV-22.  One  of  the  most  striking  features  for  this  event  was  clear 
superiority  of  detector  performance  in  the  1.5-2.  5 Hz  passband.  The  lower 
detector  output  in  the  3.  0-4,  0 Hz  passband  than  in  the  1.  5-2.  5 Hz  passband  can 
be  attributed  to  better  noise  coherency  in  the  3.  0-4.  0 Hz  passband  within  a sub- 
array and  noise  incoherency  between  the  subarrays.  Evaluating  the  detector 
performance  on  the  basis  of  the  detection  thresholds  used  at  NORSAR,  as  the 
incoming  signal  becomes  weaker,  the  detector  output  reaches  those  threshold 
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MEA^’  SURARRAY  SNR  VERSUS  DETECTOR  OUTPUT 

(WRS/248/07N  ) 


Since  the 


levels  in  the  order:  curves  D^,  C^,  A^,  D^,  C^,  and 

coherent  and  incoherent  detectors  in  the  1.5-2.  5 Hz  passband  have  11  dI3  and 
4 dB  output,  respectively,  at  the  9 dB  level  of  mean  subarray  SNR,  the  above 
order  is  not  sufficiently  significant  to  resolve  the  difference  in  performance 
of  the  detectors. 

In  general,  the  event  KAZ/282/06N  has  better  signal  similarity 
than  the  event  WRS/248/07N  because  the  former  event,  a presumed  explosion, 
yields  higher  detector  output  than  the  latter,  an  earthquake,  at  the  same  level 
of  mean  subarray  SNR.  The  presumed  explosion  showed  more  energy  in  the 
3.  0-4.  0 Hz  passband  than  the  earthquake  event  as  well  as  higher  signal  coher- 
ency. 

3.  The  Measured  Detection  Probability 

Figures  lV-24  through  IV -27  show  the  measured  detection  pro- 
bability from  the  ensemble  of  91  events  for  various  cases.  Figures  lV-24  and 
IV-25  are  for  the  updated  LTAmode,  using  (STA/LTA)  , in  the  1.  5-2.  5 

Hz  and  3.  0-4.  0 Hz  passbands,  respectively,  and  Figures  IV-26  and  lV-27  for 
the  frozen  and  averaged  LT A mode,  using  (STA)  / LTA  , in  the  same  two 
passbands. 

For  a given  decision  threshold,  the  diversity-stack  beams  show 
a higher  detection  probability,  particularly  for  the  incoherent  detectors.  How- 
ever, for  weaker  signals  the  differences  between  the  adjusted-delay  and  the 
diversity-stack  beams  are  smaller  because  the  signal  amplitudes  are  compar- 
able to  the  noise  level  and  thus  the  channel  weights  are  more  random. 

The  detection  probability  is  highly  dependent  upon  the  distril)u- 
tion  of  signal  amplitude  involved.  As  mentioned  earlier,  the  event  selection 
was  constrained  to  the  NORSAR  near -regional  earthquakes  ai*.d  the  presumed 
explosions  available  on  the  existing  NORSAR  short-period  subarray  beam  tapes. 
The  bodywavc  magnitudes  were  not  uniformly  distributed  for  ti'.-;se  data  because 
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the  36  presumed  explosions  were  of  larger  magnitude  and  generally  had 
stronger  P-waves  (see  Table  III-2).  The  resulting  probability  curves  for 
both  the  coherent  and  incoherent  detectors  are  significantly  affected  by 
these  strong  events,  particularly  for  decision  thresholds  greater  than  18  dB 
and  36  dB  levels  for  the  updated  and  frozen  LTA  modes,  respectively. 

In  brief,  tlie  diversity-stack  beams  yield  higher  detection 
probability  for  both  coherent  and  incoherent  detectors  as  generally  expected. 
The  next  subsection  presents  the  evaluation  of  these  detectors  on  the  basis 
of  their  operating  cliaracteristic  curves. 

E.  THE  MEASURED  DETECTOR  OPERATING  CHARACTERISTICS 

In  this  subsection  the  false  alarm  rate  and  detection  probability 
measurements  are  combined  together  to  form  the  detector  operating  character- 
istics or  operating  curves  for  constant  detector  output.  The  conventional 
evaluation  of  detector  performance  can  be  accomplished  by  appraising  these 
clia rac ter i s tic s.  In  the  last  two  subsections,  it  has  been  shown  that  the 
diversity-stack  beams  outperform  the  adjusted-delay  beams  for  signal 
detection  and  have  more  uniform  false  alarm  probability  distribution  than 
the  adjusted  delay  beams  for  both  coherent  and  incoherent  detectors.  Con- 
sequently, tlie  operating  curves  are  formed  for  the  diversity- stack  beams  only. 

Figures  IV-12  and  IV- 1 3 (False  Alarm  R-obability  for  Diversity - 
Stack  Beams)  and  Figures  IV- 24  to  IV-27  (Detection  Probability  for  Diver- 
sity-Stack Beams)  indicate  that  the  overlapping  range  of  tlie  actual  detector 
output  is  not  sufficiently  large  enough  to  form  the  characteristic  functions  by 
c.lirninating  the  constant  delector  outfiut  (signal-to-noise  ratio  in  dB)  because 
the  events  used  in  this  study  usually  have  strong  signals  and  the  detector  outiiut 
in  noise  is  very  low.  One  way  to  surmount  this  difficulty  is  to  extrapolate 
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the  false  alarm  probability  curves  in  Figures  IV-12  and  IV-13  to  the  higher 
thresholds  sufficiently  high  enough  to  form  the  operating  curves  with  the 
detection  probabilities  in  Figures  IV- 24  to  IV-27  • The  parameters  of 
Gaussian  distribution  functions,  mean  values  and  standard  deviations,  for 
the  observed  false  alarm  probabilities  in  Subsection  IV-C  can  be  used  for 
this  purpose.  For  the  probabilities  lees  than  10  of  normal  distribution 
functions,  the  mathematical  handbook  compiled  by  Abramowitz  and  Stegun 
(1966)  can  serve  as  a good  aid  for  approximate  computation. 
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Figure  IV-28  shows  the  operating  characteristics  for  the 

updated  LTA  mode,  using  (STA/LTA)  for  the  coherent  and  incoherent 

max 

detector  s in  the  1.  5-2.  5 Hz  and  3.  0-4.  0 Hz  passbands,  respectively.  The 

-9 

false  alarm  probabilities  have  been  extrapolated  to  the  10  levels  (17.  5 dB 

and  15.8  dB  for  the  coherent  detector  in  the  1 . 5-2.  5 Hz  and  3.  0 -4.  0 Hz 

passbands,  respectively,  and  4.7  dB  and  5.6  dB  for  the  incoherent  detector 

in  the  1.5-2.  5 Hz  and  3.  0-4.0  Hz  passbands,  respectively)  and  are  plotted 

in  logarithms  while  the  detection  probabilities  are  shown  in  linear  scale. 

For  the  same  level  of  false  alarm  probabilities,  the  detector  performance 

is  ranked  in  the  order:  incoherent  detector  operating  in  the  1.  5-2.  5 Hz  pass- 

band,  coherent  detector  in  the  1.  5-2.  5 Hz  passband,  incoherent  detector  in 

the  3.  0-4.0  Hz  passband  and  coherent  detector  in  the  3.  0-4.0  Hz  passband. 

This  fact  means  that  the  incoherent  detectors  outperform  the  coherent 

detectors  for  each  passband  and  it  is  more  pronounced  in  the  range  of  low 

false  alarm  rate  because  the  long-term-average,  LTA,  for  the  incoherent 

beam  increases  slower  than  the  coherent  beam  for  the  (STA/LTA)  mode 

max 

ol  detectors  after  the  onset  of  signal  arrivals.  For  the  weak-signal  range, 
it  still  clearly  remains  true  that  the  incoherent  detector  yields  better  per- 
formance than  the  coherent  detector  for  the  updated  LTA  mode  which  is  the 
"Operational"  mode  at  NORSAR  on-line  detection  system. 
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Figure  IV-29  shows  the  detector  operating  characteristics  for 

the  frozen  and  averaged  LTA  mode,  using  (STA)  /LTA,  for  the  diversity- 

max  ' 

stack  beams  in  the  1.5-2.  5 Hz  and  3.0-4.  0 Hz  passbands.  For  the  strong 
signal  range  and  low  false  alarm  rate  level,  the  incoherent  detector  clearly 
demonstrates  the  superior  performance  over  the  coherent  detector  in  both 
passbands.  For  the  3.  0-4.0  Hz  passband,  the  coherent  and  incoherent  detectors 
yield  about  the  same  performance  with  slight  indication  of  incoherent  detector 
being  better  than  the  coherent  detector.  However,  for  the  1.  5-2.  5 Hz  pass- 
band  , both  coherent  and  incoherent  detectors  show  the  same  performance  for 
the  present  ensemble  of  events.  Finer  resolution  of  detection  probability 
in  the  weak-event  range  might  tell  the  difference. 

In  summary,  from  the  evaluation  of  detector  operating  character- 
istics, it  is  concluded  that  the  incoherent  detector  operating  in  the  1.  5-2.  5 H/. 
passband  yields  best  performance  and  the  coherent  detector  in  the  3.  0-4,0  Hz 
passband  the  worse  among  the  four  cases  studied.  For  each  passband,  the  incoherent 
detector  shows  the  better  performance  than  the  coherent  detector.  For 
detecting  the  underground  explosions  and  near-regional  earthquakes,  the 
incoherent  detector  is  superior  to  the  coherent  detector. 
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CONCLUSIONS 


This  section  summarizes  the  major  results  and  conclusions  of 


this  study. 


A. 


FALSE  ALARM  RATE  MEASUREMENTS 


The  false  alarm  probability  density  function  from  either  the  co- 
herent or  the  incoherent  detector  output  closely  follows  a Gaus- 
sian distribution.  This  indicates  that  the  amplitude  distribution 
of  the  envelope  of  the  noise  beam  can  be  modeled  by  a lognormal 
distribution. 


For  the  coherent  envelope  detector,  the  1.5-2.  5 Hz  passband  had 
a higher  false  alarm  probability  than  the  3.  0-4.  0 Hz  passband  for 
a given  decision  threshold. 


For  the  incoherent  detector,  the  opposite  is  true.  For  a given 
decision  threshold,  the  3.  0-4.  0 Hz  passband  bad  a higher  false 
alarm  probability  than  the  1. 5-2.  5 Hz  passband. 


The  detectors  using  true  envelopes  computed  with  Hilbert- 
transform  yielded  higher  false  alarms  than  the  STA-envelope 
detectors  because  the  former  did  not  incorporate  any  smoothing 
through  integration. 


In  order  to  achieve  the  false  alarm  rates  less  than  one/per 
beam/per  day,  the  decision  threshold  for  detection  must  be  lar- 
ger than  9 dB  for  the  coherent  detector  in  both  passbands.  For 
the  incoherent  detector,  the  thresholds  are  3 dB  in  the  1.5-2.  5 Hz 
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passband  and  4.  5 dB  in  the  3.  0-4.  0 Hz  passband. 


13.  DETECTION  PERFORMANCE 

In  this  subsection,  the  concluding  summary  on  detection  is 
l)ased  only  on  the  detector  output  for  the  presumed  explosions  and  near- 
regional  earthquakes  without  taking  into  account  the  detector's  false  alarm 

rates.  Overall  conclusions  on  relative  detector  performance  are  in  Subsection 
V-C. 

For  the  operational  mode  detectors  using  a running  STA/LTA 
ratio,  the  coherent  detector  output  was  2-3  dB  higher  than  the 
incoherent  detector  in  both  passbands  for  the  presumed- explosion 
ensemble.  For  the  earthquake  ensemble,  it  was  5-6  dB  higher. 

• For  the  frozen- LTA-mode  detectors  using  (STA)  /LTA  the 

max  ’ 

coherent  detector  yielded  a 10-11  dB  higlicr  output  than  the  in- 
coherent detector  in  the  1,5-2.  5 Hz  passiiand  and  7-8  dB  in  the 
3.0-4.  0 Hz  passband. 

• The  average  detector  output  for  dive rsity- stack  beams  using 

Hilbert  transform  pairs  was  7.  3 dB  and  4.  1 dB  higher  than  the 

STA -envelope  average  for  the  (STA/LTA)  and  (STA)  /iTtA 

iiiax  max 

coherent  detectors,  respectively,  and  5.  1 dI3  and  2.  1 dB  higher 
for  the  (STA/LTA)  and  (STA)  /iTTA  incoherent  detectors. 

lllctX  llltlX  * 

respi:cti vely , for  four  events. 

• The  simulation  study  wlurre  the  scaled  signal  amplitudes  were 
buried  in  noise  suggested  I hat  the  coherent  detector  operating 
in  the  1.5-2.  5 Hz  passband  yiebbul  th<>  best  iletection  among  the 
various  ilete<l()rs  studied.  However,  for  the  3.  0-4.  0 Ilz  passband 
th(*  i neohe  ri'ut  detector  yielded  Ix'ttei'  pi'rformance  than  the 
herent  fli‘tc«  tor 
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• A detector  where  the  LTA  is  delayed  relative  to  the  STA  or 
E(t)  envelope  can  reduce  signal  contamination  in  the  LTA  and 
thereby  provide  a better  signal-to-noise  ratio  measurement, 

• The  diversity-stack  beams  had  a higher  detection  probability 
than  the  adjusted-delay  beams  for  both  coherent  and  incoherent 
detectors. 


C.  DETECTOR  OPERATING  CHARACTERISTICS 

Final  conclusion  on  detector  performance  is  sumirarized  as  fol- 


lows: 


The  incoherent  detector  yields  better  p(*rformance  than  the  co- 
herent detector  in  both  the  1.  5-2.  5 Hz  and  the  3.  0-4.  0 Hz  pass- 
bands. 


• Both  coherent  and  incoherent  detectors  yield  better  perfoi  nance 
in  the  1.  5-2.  5 Hz  passband  than  in  the  3.  0-4.  0 Hz  passband. 

• The  envelope  detector  using  Hilbert  transformation  operating  at 
a given  level  produced  a larger  false  alarm  rate  but  an  even  bet- 
ter detection  performance  than  the  STA  type  of  detector  at  the 
equivalent  level.  Thus  the  achial  operating  characteristic  of  the 
Hilbert  transform  envelope  detector  was  superior  to  the  STA 
envelope  detector. 

In  conclusion,  for  detecting  underground  explosions  and  near- 
regional  seismic  signals  at  NORSAR,  the  incoherent  beamforming  envelope 
detector  is  superior  to  the  coherent  beamforming  envelope  detector,  as  expected. 
For  teleseismic  events,  however,  it  may  not  be  expected  that  the  incoherent 
detector  would  still  yield  the  better  performance.  Hence  it  is  worthwhile  to  use 
the  incoherent  detector  as  a supplement  to  the  coherent  detector  in  order  to 
maintain  good  detection  performance  for  events  at  all  distances. 
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